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ABSTRACT 

We present Chandra, XMM-Newton and RXTE observations of 1WGA J1713. 4-3949, a compact source 
at the center of the galactic supernova remnant (SNR) G347.3-0.5. The X-ray spectrum of the source 
is well-fitted by the sum of a blackbody component with a temperature of ~ 0.4 keV plus a power law 
component with photon index ~ 4. We found no pulsations down to 4% in the 0.01-0.16 Hz range and 
down to 25% in the 0.01-128 Hz range. This source resembles other compact central objects (CCOs) in 
SNRs, and we suggest that 1WGA J1713.4-3949 is the associated neutron star for G347.3-0.5. We also 
measured the properties of the adjacent radio pulsar PSR J1713-3945 with a 392 ms period and show 
that it is not associated with 1WGA J1713. 4-3949 nor, most probably, with SNR G347.3-0.5 as well. 

Subject headings: ISM: individual (G347.3-0.5) — pulsars: individual (J1713-3945) — stars: neutron 
(1WGA J1713. 4-3949) — supernova remnants — X-rays: stars 



(Lazendic et al. 2003). In this paper we present a detailed 
spectral and timing analysis of 1WGA J1713. 4-3949 using 
our Chandra data, as well as archival data from the XMM- 
Newton and RXTE satellites. We also present properties 
of the radio pulsar PSR J1713-3945, which imply it is not 
related to 1WGA J1713.4-3949 or SNR G347.3-0.5. 

2. OBSERVATIONS 

Two fields towards G347.3— 0.5 were observed on 2000 
July 25 with the Advanced CCD Imaging Spectrometer 
(ACIS) detector on board the Chandra X-ray Observatory 
(ObsID 736 and 737). The ACIS-I detector for ObsID 
736 was positioned at the bright northwestern rim of the 
SNR and the roll angle was such that 1WGA J1713.4- 
3949 fell on the SI chip (see Figure 1). Data were taken in 
timed exposure mode with an integration time of 3.2 s per 
frame. The fully processed Level 2 data were reduced us- 
ing standard threads in the CIA0 software package version 
2.2.1. The effective exposure time after data processing 
was 29.6 ks. The exposure-corrected Chandra image of 
the SNR and the compact source 1WGA J1713.4-3949 is 
shown in Figure 1, with the ROSAT contours overlaid. 
The image was binned in 4'.'6 x 4'.' 6 pixels and smoothed 
with a Gaussian filter with a FWHM of 2". The image 
of the compact source is elongated because its location is 
well off-axis (~ 23') from the aim point. 

We obtained archival XMM data of 1WGA J1713.4- 
3949 carried out on 2001 March 2 with a 15 ks exposure 
(ObsID 0093670501; see also Cassam-Chenai et al. 2003). 
The European Photon Imaging Camera (EPIC)-pn cam- 
era was operated in Extended Full Frame mode with time 
resolution of 199 ms, while the two EPIC-MOS cameras 
were operated in Full Frame mode with a time resolution 
of 2.6 s. The EPIC data reduction was performed us- 
ing the SAS software package version 5.4.1. The event files 
used for analysis were created from observational data files 
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1. INTRODUCTION 

For many years it has been well known that pulsars — 
compact remains of supernova explosions — are rapidly 
rotating, highly magnetized neutron stars (NSs). How- 
ever, sources discovered in recent years have shown that 
NSs also appear to come in other varieties: even more 
highly magnetized objects with slow spins, such as anoma- 
lous X-ray pulsars (AXPs) and soft gamma-ray repeaters 
(SGRs) (e.g., Mereghetti 1998; Thompson 2000), and more 
poorly understood objects near the centers of SNRs with 
detectable X-ray flux, but no optical/radio counterparts or 
signs of rotation, so-called compact central objects (CCOs; 
Pavlov et al. 2002). Their spectra are very soft, with 
blackbody temperatures ~ 0.4 keV, often requiring ad- 
ditional power law component with photon indices of ~ 4 
(Pavlov et al. 2002). In this article we address the nature 
of 1WGA J1713. 4-3949, a compact source located at the 
center of the SNR G347.3-0.5 (RX J1713. 7-3946). Early 
ROSAT observations identified two point sources within 
the boundaries of the SNR shell (Pfcffcrmann & Aschen- 
bach 1996, ; see also Figure 1). The source further away 
from the center of the remnant, 1WGA J1714.4— 3945, was 
determined to be of stellar origin (Pfeffermann & Aschen- 
bach 1996). For the other source, 1WGA J1713.4-3949, 
located at the geometrical center of the SNR, no optical 
counterpart has been found within 10" of the ROSAT po- 
sition (Slane et al. 1999), and the source has correspond- 
ingly been suggested to be a candidate for an associated 
neutron star. While ASCA observations were not able to 
provide useful limits to the pulsations in the X-ray band 
(Slane et al. 1999), radio pulsations were detected with a 
392 ms period within a 7' radius region towards this source 
(PSR J1713-3945 Crawford et al. 2002). 

1WGA J1713. 4-3949 has been observed with the Chan- 
dra X-ray Observatory as a part of a project to study 
nonthermal radio and X-ray emission from the SNR shell 
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(ODFs) using the SAS tasks epchain and emchain. The 
net exposure time after filtering event files for good time 
intervals was 9.2 ks in the pn camera and 14.1 ks in each 
MOS camera. 

We used archival data of SNR G347.3-0.5 (Pannuti ct al. 
2003) from the Proportional Counter Array (PCA) aboard 
the RXTE satellite, which is sensitive to X-ray photons in 
the energy band 2-60 keV and has a field of view with 
radius « 1°. While the PCA has no imaging capability, 
it has a large effective area (6000 cm 2 ) which makes it 
suitable for searching for faint X-ray pulsations. Obser- 
vations of G347.3-0.5 were made in 7 separate intervals 
between 1999 June 14-24, with a total exposure time of 
100 ks. Analysis of the RXTE data was carried out using 
LHEAS0FT software version 5.2. 

A 392-ms radio pulsar, PSR J1713— 3945, was discov- 
ered in a 20 cm targeted search for pulsations in SNR 
G347.3-0.5 using the Parkes telescope and multibeam re- 
ceiver. The details of the search and discovery are outlined 
in Crawford ct al. (2002). The pulsar was subsequently 
timed with the multibeam receiver at Parkes to obtain as- 
trometric and spin parameters and to determine whether 
the pulsar was associated with 1WGA J1713. 4-3949. 

3. ANALYSIS AND RESULTS 

3.1. Spectral Analysis 

For the spectral analysis of the Chandra data, we ex- 
tracted the source counts from a 37" x 61" elliptical re- 
gion around the compact source, and for the background 
region we used a 2' x 3' elliptical region north of the 
source. We obtained 8454 ± 103 background-subtracted 
counts for the source and the spectrum was grouped to 
contain at least 25 counts per bin. Residuals from fit- 
ting standard continuum models to the source spectrum 
(blackbody, bremsstrahlung and power law) revealed an 
apparent change in the gain calibration on chip SI that 
results in a spurious absorption feature around 2 keV. We 
extracted data from the on-board calibration source taken 
before and after our observation, and these show that there 
is indeed a gain shift for this CCD of ~ 3.5%. However, we 
found that the fit parameters did not change significantly 
when the gain was frozen at this 3.5% value or when it 
was left as a free parameter. We thus allowed gain to vary 
in our models (see Table 1). 

For spectral analysis of the XMM data, we obtained 
9600±217 background-subtracted counts from the pn cam- 
era, and 5451 ± 74 (5274 ± 73) in the MOS1 (MOS2) cam- 
era. In the pn camera, the counts were extracted from a 
27" radius circular region around the source, and for the 
background extraction we used a 3'4 x 8'3 box north from 
the source position inside the same CCD chip. Similarly, 
in the MOS1 (MOS2) camera we used 22" (19") radius 
circular region for the source and 3f4 (2.'2) radius circular 
region on the same chip for the background. The spectra 
were also grouped to contain at least 25 counts per bin. 
We performed a joint spectral analysis of the three EPIC 
data sets, allowing only the relative normalization to vary. 
The results obtained were consistent with those obtained 
from the Chandra spectrum. We therefore proceeded with 
a joint spectral analysis of Chandra and XMM spectra. 

The results from fitting a few typical neutron star mod- 

5 http://pulsar.princeton.edu/tempo 



els are summarized in Table 1 and a fit example is given 
in Figure 2. Fits with the bremsstrahlung and black- 
body models gave equally acceptable results, while the 
power law model gave a somewhat higher reduced x 2 - The 
steep photon index of T ~ 4 excludes the possibility that 
1WGA J1713. 4-3949 is a background active galactic nu- 
cleus (AGN), which typically have photon indices of 1.2- 
2.2 (Turner & Pounds 1986). The best-fit column density 
(Njj) in the power law model of ~ llxl0 21 cm~ 2 is also 
somewhat higher than the value of 8 x 10 21 cm~ 2 found to- 
wards the SNR (Lazendic et al. 2003). On the other hand, 
a blackbody model yields Nh ~ 4xl0 21 cm~ 2 , which is a 
lower value than found towards the SNR. The blackbody 
fit is improved by adding a power law component, which 
also results in an Nh value consistent with that found 
towards the SNR. In this composite spectrum the unab- 
sorbed source luminosity is Lx(0.5 — 10.0 keV) ~6x 10 34 
erg s _1 , adopting a distance of ~ 6 kpc to the SNR (Slane 
et al. 1999). The blackbody normalization yields a rela- 
tively small emitting area with radius R ~ 2.4 £>6fcpc km, 
which is significantly smaller than the canonical value for 
the NS radius of 10 km. This could imply that most of the 
thermal emission originates from the polar caps on the NS 
surface, while the thermal component from the rest of the 
NS surface is too faint to be detected. 

The small inferred radius may be due, at least partly, 
to oversimplification of the thermal emission from the NS 
by the blackbody model. More appropriate models that 
take into account the NS atmosphere (NSA) yield lower 
NS temperatures, which in turn correspond to larger NS 
radii (e.g., Zavlin, Pavlov, & Shibanov 1996). To investi- 
gate this, we fit our spectra with two NSA models. The 
NSA model by Gansicke, Braje, & Romani (2002) is calcu- 
lated for weakly magnetic NS (B < 10 10 G) and considers 
hydrogen, solar or iron compositions of the atmosphere. 
The NSA model by Pavlov et al. (1995) assumes the stan- 
dard magnetic field strength (B = 10 12 G) and hydrogen 
atmosphere. Both models assume a standard NS mass 
of 1.4M Q . Both hydrogen models gave acceptable residu- 
als and yielded emitting region sizes close to 10 km, with 
corresponding reductions in the NS temperature. For ei- 
ther model, however, the derived temperature is above the 
expected surface temperature for a cooling neutron star 
more than several hundred years old. Model atmospheres 
with solar abundance or iron composition gave unaccept- 
able fits predicting broad absorption features not present 
in the spectra. 

3.2. Timing Analysis 

We used a Fast-Fourier Transform (FFT) to search the 
X-ray data for pulsations from 1WGA J1713.4-3949, but 
found no pulsed signal; we list upper limits on the pulsed 
fraction and the search parameters in Table 2. The up- 
per limit from RXTE data is least constraining because 
1WGA J1713.4-3949 only accounts for 0.6% of the total 
flux in the RXTE's large field of view. 

To derive the timing properties of the radio pulsar PSR 
J1713-3945 the TEMPO timing package 5 was used on a total 
of 27 times-of-arrival (TOAs) with the JPL DE200 plane- 
tary ephemeris (Standish 1990). The TEMPO package con- 
verts each TOA to the solar system barycenter then refines 
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the timing parameters by minimizing residuals between 
observed and model TOAs over the observation span. The 
resulting timing parameters are presented in Table 3. The 
pulsar's dispersion measure (DM) of 337 pc cm~ 3 yields 
an estimated distance of 4.3 kpc using the revised DM- 
distance model of Cordes & Lazio (2002). 

4. DISCUSSION 

The position of the 1WGA J1713.4-3949 determined 
from the XMM data is R.A (J2000) 17 h 13 m 28 s .4, Decl. 
(J2000) -39°49'54'/5 (with an uncertainty of - 6"). Us- 
ing our 20 cm radio data of SNR G347.3-0.5 obtained 
with Australia Telescope Compact Array (Lazendic et al. 
2003), we found no radio counterpart for 1WGA J1713.4- 
3949, and derive a 5c upper limit for a point source flux 
of 3 mJy. This is similar to the limit derived towards 
the CCO in Puppis A (Gaensler, Bock, & Stappers 2000). 
There is no apparent optical emission at the location of 
1WGA J1713.4-3949 in the Digitized Sky Survey 6 image, 
which has an upper limit of V > 17 mag. 

While the distance to the radio pulsar PSR J1713-3945 
is broadly consistent with the SNR distance of 6 ± 1 kpc 
estimated by Slane et al. (1999), the pulsar is likely to be 
quite old (r c ~ 1.1 Myr, compared to t < 40 kyr for SNR 
G347.3-0.5) and does not have sufficient spin-down lumi- 
nosity (E ~ 3.7x 10 33 erg s _1 ) to power the observed X-ray 
flux from 1WGA J1713.4-3949 (L x ~ 6 x 10 34 erg s^ 1 ). 
The newly determined radio timing position of the pul- 
sar is ~ 4' in declination away from the position of 
1WGA J1713.4-3949 (see Figure 1). PSR J1713-3945 
is therefore spatially coincident with SNR G347.3-0.5 by 
chance, and there is likely no physical association between 
the two systems. We detected no X-ray source at the posi- 
tion of PSR J1713-3945 in Chandra or XMM data. Adopt- 
ing Njj ~5x 10 21 cm~ 3 (by scaling the value for the SNR 
to the 4.3 kpc distance of this pulsar), we derive an up- 
per limit on the luminosity in the 0.5-10.0 keV band of 
2 x 10 31 erg s _1 for nonthermal emission with a photon 
index ~ 1.5, and an upper limit on the temperature of 
82 eV for blackbody emission from a 10-km radius NS. 

The lack of radio and optical counterparts for 
1WGA J1713. 4-3949, the absence of X-ray pulsations 
with the current sensitivity and time resolution, the two- 
component spectrum and its associated luminosity are 
properties consistent with CCOs (Pavlov et al. 2002). 

If this CCO is an accretion powered NS, we can consider 
two scenarios: accretion from a fallback disk or a low- 
mass companion. The fallback disk has been suggested to 
form from residual material after the SNR explosion (e.g., 
Chevalier 1989) and has been used to explain and unify the 
properties of various types of NSs (Alpar 2001). However, 



limits derived from optical and infrared observations seem 
to disfavor this model, at least for some types of NSs (e.g., 
Kaplan et al. 2002). X-ray monitoring, in conjunction with 
optical and infrared limits, can be used to establish if the 
NS is variable and thus likely to be accreting from a low- 
mass companion. Indeed, variability and a period of 6.4 h 
has been reported for one of the CCOs, J1617-5102 in 
SNR RCW 103, on which basis this CCO is claimed to 
be a binary (Sanwal et al. 2002). Comparing the unab- 
sorbed X-ray flux of 1WGA J1713.4-3949 obtained from 
Chandra and XMM data for a blackbody fit with that of 
ASCA data (F&(0.5 - 10.0) = 5.3xl0~ 12 ergem"^" 1 ; 
Slane et al. 1999), we find no long-term (~ 3 yr) variabil- 
ity in 1WGA J1713.4-3949. 

Interpretations of the emission from CCOs as thermal 
radiation from the surface of a NS run into problems due 
to high surface temperatures and small emitting areas de- 
rived from the blackbody model fits. The temperature 
we derived from the blackbody model of 0.4 keV is much 
higher that the temperature predicted by the standard 
cooling model (Page & Sarmiento 1996), while the tem- 
perature derived from the NSA models (0.2 keV) would 
imply an NS age of 100-1000 yr, which is less than the 
estimated 20-40 kyr age of the SNR (Slane et al. 1999). 

In the case that CCOs are rotation powered NSs, the 
non-detection of radio pulsations can be explained by 
unfavorable beaming orientation. One of the objects, 
IE 1207.4-5209 in SNR G296.5+10.0, that was originally 
identified as a CCO, has been found to have a period of 
424 ms with faint pulsed fraction (7.6%; Zavlin et al. 2000). 
Our limits on pulsed fraction from 1WGA J1713.4-3949 
are consistent with those derived towards other CCOs (7- 
15%; e.g., Pavlov et al. 2002). Our limits indicate that 
1WGA J1713.4-3949 is probably not an AXP, which are 
generally found to have spin periods between 6-12 s and 
pulse fractions of 10-70% (e.g., Mereghetti 1998). 

In summary, the nature of the compact central object 
in SNR G347.3-0.5 remains uncertain. Future observa- 
tions are planned to help determine the exact nature of 
this source. Pointed Chandra observations arc needed to 
provide a more accurate position and to enable a search for 
an optical counterpart, while deeper timing observations 
can search for periods substantially shorter than available 
from the current data. 



We thank Peter Woods for useful discussions and Slava 
Zavlin for providing the files for the magnetic neutron 
star atmosphere. This work was supported in part by 
NASA contract NAS8-39073 (POS), grant GO0-1123X 
(JSL) and grant GO2-3080B (JPH). 



REFERENCES 



Alpar, M. A. 2001, ApJ, 554, 1245 

Cassam-Chcnai, G., Dccourchelle, A., Ballet, J., Sauvageot, J.-L., 

Dubner, G., 2003, A&A, in preparation 
Chevalier, R. A. 1989, ApJ, 346, 847 

Cordes, J. M. & Lazio, T. J. M. 2002, preprint, astro-ph/0207156 
Crawford, F., Pivovaroff, M. J., Kaspi, V. M., & Manchester, R. N. 

2002, ASP Conf. Ser. 271: Neutron Stars in Supernova Remnants, 

cds. P. O. Slane & B. M. Gaenlscr, 37 



Gaensler, B. M., Bock, D. C.-J., & Stappers, B. W. 2000, ApJ, 537, 
L35 

Gansickc, B. T., Braje, T. M., & Romani, R. W. 2002, A&A, 386, 
1001 

Gotthclf, E. V., Pctrc, R., & Vasisht, G. 1999, ApJ, 514, L107 
Kaplan, D. L., Kulkarni, S. R., Frail, D. A., & van Kerkwijk, M. H. 

2002, ApJ, 566, 378 
Lazendic, J., Slane, P., Gaensler, B., Plucinsky, P., Hughes, J. & 

Reynolds, S., 2003, ApJ, submitted 



DSS was produced at the Space Telescope Science Institute under U.S. Government grant NAG W-2166. 



4 



Mereghctti, S. 1998, Memoric dclla Societa Astronomica Italiana, 69, 
819 

Pannuti, T., Allen, G., Houck, J., Sturner, S. 2003, ApJ, 593, in press 
Page, D. & Sarmiento, A. 1996, ApJ, 473, 1067 

Pavlov, G. G., Shibanov, Y. A., Zavlin, V. E., & Meyer, R. D. 1995, 

in: The Lives of the Neutron Stars, eds. M. A. Alpar, U. Kiziloglu, 

& J. van Paradijs (Kluwer: Dordrecht) p. 71. 
Pavlov, G. G., Sanwal, D., Garmire, G. P., & Zavlin, V. E. 2002, 

ASP Conf. Ser. 271: Neutron Stars in Supernova Remnants, eds. 

P. O. Slane & B. M. Gaenlser, 247 
Pfeffermann, E. & Aschcnbach, B. 1996, Rontgcnstrahlung from the 

Universe, eds. Zimmermann, H.U., Trumper, J., and Yorkc, H., 

MPE Report 263, 267 



Sanwal, D., Garmire, G. P., Garmire, A., Pavlov, G. G., & Mignani, 

R. 2002, American Astronomical Society Meeting, 200, #72.01 
Standish, E. M. 1990, A&A, 233, 252 

Slane, P., Gacnsler, B. M., Dame, T. M., Hughes, J. P., Plucinsky, 

P. P., & Green, A. 1999, ApJ, 525, 35 
Thompson, C. 2000, ASP Conf. Ser. 202: IAU Colloq. 177: Pulsar 

Astronomy - 2000 and Beyond, 669 
Turner, T. J. & Pounds, K. A. 1986, MNRAS, 240, 833 
Vaughan, B. A. et al. 1994, ApJ, 435, 362 

Zavlin, V. E., Pavlov, G. G., & Shibanov, Y. A. 1996, A&A, 315, 
141 

Zavlin, V. E., Pavlov, G. G., Sanwal, D., & Trumper, J. 2000, ApJ, 
540, L25 



5 



Table 1 

Parameters of the joint spectral fits to the Chandra and XMM spectra of the compact source 

1WGA J1713. 4-3949 with different models. 



Model 


N H 


kT 


r 


r? a 


Xljdoi Gain b 


po 

r X(0. 5-10.0) 




(10 21 cm- 2 ) 


(keV) 




(km) 


(%) 


( ergcm~ 2 s _1 ) 



BREMSS 




Q2+ - 02 






1.07/675 


c 1+0.5 
°- 1 -0.2 


9.4xl0~ 12 




POW 






4 2+ ' 1 




1.44/675 


4.6l8: 8 5 


3.2X10" 11 




BB 


o fi +0.1 


n 4n+ 001 




2 5+ ' 1 


1.25/675 


6.018'j 


4.5xl0~ 12 




POW+BB 


o n+0-9 

°- U -0.6 


°- 38 -0^02 


o q+0-3 

O.»_ 2 


2 4+0.6 


1.05/670 


7 1+0.9 


1.4(0.3)xl0" 


n 


POW+NSA(G) d 


7 o+0.6 
' •' 3 -0.4 


o.i8l8:S 


o 7+O.8 
°- ' -0.3 


13.0+1° 


1.07/670 


c c+0.6 
°-°-0.5 


l.l(0.4)xl0" 


n 


POW+NSA(Z)° 


°' 1 -0.9 


o.26lS:g 


o n+0.2 

^-0.5 




1.22/667 


c 1+0.5 
°- 1 -0.5 


1.5(0.3)xl0" 


li 



a radius of the emitting region at infinity assuming the source distance of 6 kpc 
b a linear gain shift included for fitting the Chandra data 

c unabsorbed flux — for the two component models the flux of the second component alone is given in 
parentheses 

d NSA model by Gansicke, Braje, & Romani (2002) for weak magnetic field of < 10 10 G 
°NSA model by Pavlov et al. (1995) for standard magnetic field of 10 12 G 



Table 2 

Parameters and results of pulsation search from 1WGA J1713. 4-3949 using FFT. 



Detector 


Time 


Frequency range 


f p limit a 




resolution 


searched 




ACIS-S1 


3.2 s 


0.01-0.16 Hz 


12% 


EPIC-MOS b 


2.6 s 


0.01-0.19 Hz 


12% 


EPIC-pn 


199 ms 


0.01-2.5 Hz 


4% 


PCA C 


4 ms 


0.01-128 Hz 


25% 



fractional pulsation upper limit to 99% confidence level 
calculated according to Vaughan et al. (1994) 

b from combined MOS1 and MOS2 data 

c calculated from FFT of two separate intervals of <~ 3 
days of data with gaps 



Table 3 

Properties of the radio pulsar PSR J1713-3945 detected towards G347.3-0.5. 



Parameter Valuc a 



R.A. (J2000) 17 h 13 m 14 s 23(3) 

Decl. (J2000) -39°45'47'.'5(14) 

Period, P (ms) 392.451769234(8) 

Period derivative, P (xl0~ 15 ) 5.7380(5) 

Dispersion Measure, DM (pc cm" 3 ) 337(3) 

Epoch of period (MJD) 52165.0 

rms residual (ms) 1.1 

Number of TOAs 27 

Timing span (Days) 1200 

Characteristic age!? t c (Myr) 1.08 

Surface magnetic field? B (xlO 12 G) 1.52 

Spin-down luminosity? E (erg s" 1 ) 3.74xl0 33 



a Values in parentheses represent the la uncertainty in the 
least significant digit quoted. 

b T c EE P/2P 

C B ee 3.2 x 10 19 (PP) 1 / 2 G, where P is in seconds 
d E ee 4it 2 IP/P 3 , with I = 10 45 g cm 2 assumed 
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RA (J2000) 

Fig. 1. — The 0.5-10.0 kcV band Chandra greyscale image of the two fields observed towards the SNR G347.3— 0.5 (for more details see 
Lazendic et al. 2003). To show the whole extent of the remnant, the Chandra image is overlaid with the ROSAT PSPC contours (Slane et 
al. 1999). Two point sources detected with ROSA T observations are labeled, as well as the location of the radio pulsar PSR J1713— 3945. 
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Fig. 2. — ACIS and EPIC spectra from the compact source 1WGA J1713. 4-3949 and residuals for blackbody+power law model. The fit 
parameters are listed in Table 1. 



